Activated escape over oscillating barriers: The case of many dimensions 
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We present a novel path-integral method for the determination of time-dependent and time- 
averaged reaction rates in multidimensional, periodically driven escape problems at weak thermal 
noise. The so obtained general expressions are evaluated explicitly for the situation of a sinusoidally 
driven, damped particle with inertia moving in a metastable, piecewise parabolic potential. A 
comparison with data from Monte-Carlo simulations yields a very good agreement with our analytic 
results over a wide parameter range. 

PACS numbers: 05.40.-a, 82.20.Mj, 82.20.Pm 

I. INTRODUCTION 

Thermally activated escape problems in the presence of an explicit time-dependent driving are at the root of 
many timely transport processes. Typical examples comprise the control of chemical reactions with tailored laser 
pulses [0, H , ion transport through voltage-gated channels |^ , the pumping and shuttling of particles in Brownian 
environments p], ^ 0, ^ or the amplification of weak information-carrying signals via the phenomenon of 

Stochastic Resonance 
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I , to name only a few. In the absence of such a time-dependent driving and for the 

case of weak thermal noise the escape time is governed - as commonly known - by an exponentially leading Arrhenius 
Pioneered by Kramers and extended to arbitrary dimensions in the works mh, nS, EOl 



factor |T|, |T6|. 

this scheme, however, meets formidable difficulties under far from thermal equilibrium conditions. This is so because 
of an extremely complex interplay between the global properties of the metastable potential and the nonlinear noisy 
dynamics g2|, |2^, The subject of our present paper is one of the simplest and experimentally most natural 
such non-equilibrium descendants of Kramers' original escape problem ||l7|| , namely the thermally activated escape of a 
Brownian particle over a potential barrier in the presence of periodic driving which modulates both the corresponding 
potential well region and the activation barrier. While most previous attempts have been restricted to weak |2^, |2^, p7| , 
slow H, H, or fast [p| 



30| driving, we have addressed in recent analytical explorations |31, B2l by means of 



a path-integral technique the most challenging intermediate regime of moderately strong and moderately fast driving 
for a one-dimensional, overdamped escape problem. Closely related to our recent works are the subsequent appealing 
efforts in Ref . [3^ , wherein on uses instead the method of singular perturbation theory in the weak noise limit — 
a so termed WKB approximation — which, however, also has been restricted to cover only the overdamped, one- 
dimensional case. The scheme in Ref. |33 yields results which are consistent with the findings in Ref. |31, 3^; but it 
remains on a more formal and implicit level compared to ours. 

The objective of this study is to extend these recent works of time-dependent rate theory put forward in Refs. ]3l| , 
|3^ , ^ to the case of many dimension of the underlying stochastic process. In particular, we shall consider the generic 
case of a driven inertial Brownian motion dynamics. Furthermore, we shall derive for weak noise asymptotically exact 
results in analytically closed form for the escape of an inertial, sinusoidally driven Brownian particle in a metastable, 
piecewise parabolic potential in the regime of moderate forcing strengths and forcing frequencies. A comparison with 
data from Monte-Carlo simulations yields very good agreement over a wide parameter range. 

We anticipate here that the following presentation is on purpose kept rather concise. While we are confident that 
the general idea of our approach remains more transparent in this way, it is nevertheless advisable to consult the 
detailed discussion of the one-dimensional, overdamped (i.e. no inertial) dynamics in Ref. for a more thorough 
and in-depth understanding. It must be pointed out, however, that our present generalization to many dimensions in 
addition requires several conceptually new steps and features as compared to the one-dimensional overdamped case. 

In Sect. 0, we describe the model under investigation and define both the instantaneous and the time-averaged 
escape rate. The path-integral formalism and its evaluation for weak-noise is described in Sect. III. Its application 
to the escape problem is demonstrated in Sect. [V. Our general findings are then evaluated and discussed explicitly 
for the case of the driven Kramers problem with a piecewise parabolic potential in Sect. 0. Concluding remarks are 
given in Sect. VI. 
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FIG. 1: Sketch of a typical metastable potential V{x), namely the piecewise parabolic potential ( |7^ ) with AV — 0.9, muj^ — 0.6, 
and muj'^ — —0.3 in arbitrary, dimensionless units. 



II. THE ESCAPE PROBLEM 
A. The general Model 

The starting-point of our investigations is the following model for the d-dimensional Brownian motion of a particle 
with coordinates x(t) in a time-dependent force field F(x, i): 

x(i) =F(x(t),t) +\/2^B|(i) . (1) 

Here, bold quantities denote d-dimensional vectors, while d x d-matrices are represented by sans-serif fonts. The 
d-dimensional Gaussian uncorrelated, white noise ^{t) is defined by the relations (i, j = 1, 2, . . . , d) 

m)) = o, {mi,it')) = s,,sit-t') , (2) 

where Sij is the Kronecker delta and 6(t) Dirac's (5-distribution. The additive coupling of the noise to the stochastic 
dynamics is determined by the matrix B, which we assume to be non-singular. 

We restrict ourselves to the case of a time-periodic driving, i.e. F(x, i) = F(x, i + T) for a certain period T. 
Furthermore, we assume that the deterministic dynamics, i.e. Eq. for e = 0, possesses exactly one stable periodic 
orbit Xs(t) — Xs{t + T) with a time-dependent domain of attraction A{t). All other deterministic orbits, which start 
outside of A{t) and its boundary S{t) :— dA{t) are assumed to diverge in the long time limit. The boundary S{t) then 
acts as a separatrix between these two kinds of deterministic solutions. Moreover, we require that there is exactly 
one unstable periodic orbit Xtj(t) = x„(t -f T), which then moves with the separatrix, i.e. x„(i) S S{t) for all times t. 
Finally, all other deterministic solutions starting in a neighborhood of the unstable periodic orbit on the separatrix 
are assumed to be bounded. All these assumptions can usually be taken for granted in typical periodically driven 
escape problems of practical interest. 

A particularly prominent example is the Kramers problem , namely the escape of a particle with mass m out of 
the bottom well of a static potential as cartooned in Fig. ^. Under the additional infiuence of an additive sinusoidal 
driving this corresponds to the stochastic dynamics 



-V'{x) + A sin(f7 1) + y%k^£_{t) 



(3) 
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Here, rj is the viscous friction coefficient, fc^ Boltzmann's constant, T the temperature and = 2ir/T the angular 
frequency of the driving. We can rewrite the second order equation (H) in the form of Eq. , if we identify 

x:.0, F(M:^(x^,„V,„), B ^ ?) . .:-..T. (4) 

where we have introduced the one-dimensional force field F[x,t) :— —V'{x) + Asin{flt) and the frequency 7 := r]/m. 
Furthermore, we have added an auxiliary noise source of a strength proportional to S to ensure that B is non-singular. 
Eventually, we shall consider the limit 5 ^ 0, in which Eq. reduces to the original dynamics (^). Figure || depicts 
the stable and unstable periodic orbits, the domain of attraction A{t), and the separatrix S{t) for a representative 
metastable potential. 
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FIG. 2: Phase space diagram for a typical metastable potential with additive sinusoidal driving (j^, namely the piecewise 
parabolic potential ( [75[ ) with Xs = io'^ — —1, Xu = u>s = m — AV — Q — 1, rj — 0.5, A — 0.2. Solid line: Stable (left) and 
unstable (right) periodic orbit. Dashed lines: Separatrix at different times ti = 0, t2 — T/3, and ts = 2T/3. Grey area: 
Domain of attraction A{ti) of stable periodic orbit Xs{t) at time ti. 

A fully equivalent description of the stochastic dynamics (|^) is provided by the Fokker- Planck equation 1 34 for the 
probability distribution p(x, i) of an ensemble of particles x(t), 

-p(x,t)+V-j(x,i) =0 , (5) 

where the probability current density is given by 

j(x,i) :=F(x,i)p(x,t)-eDVp(x,i) (6) 
with the positive semidefinite diffusion matrix D :— BB^ . 



B. Escape rates 

If we now consider our stochastic dynamics (0) for finite noise-strengths e, it is well-known that a particle x{t) 
starting in A{t) will be able to escape out of this basin. To quantify this escape of particles, we first introduce the 
population in A{t): 

PA{t) j d'^xp(x,t) . (7) 

A(t) 
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A natural definition of an "instantaneous rate" r{t) for such escape events is then given by the relative decrease of 



this population per unit time |31, ^ 

m := -PA{t)lPA{t) . (8) 

We remark that — except for transients at early times — this quantity is independent of the initial conditions at time 
t = tf), provided we start with a distribution that is concentrated within the basin A{t). We thus make the convenient 
choice p(x, io) — (5(x — Xs(to))- Furthermore, on the time scale of these transients, T{t) will approach a time-periodic 
limit, permitting a meaningful definition of the time-averaged rate 

r:=- dt'nt'). (9) 

In the sequel of the present paper we will deduce analytic expressions for both the instantaneous and the time-averaged 
rate in the limit of small noise-strength e. 

To make progress in this direction, we first note that for weak noise e the typical time scale of the escape events 
1/r is well separated from the time scale of the just mentioned transients iQ. We can thus approximate PA{t) by 
its initial value -PA(io) = 1 in the denominator of Eq. (H). Using the Fokker-Planck equation (||) and the divergence 
theorem, the instantaneous rate reads 

T{t)^ yd'^-in.[j(x,t)-xp(x,i)] . (10) 

s(t) 

Here, the integration is over the entire time-dependent separatrix S{t) and n denotes their outer normal vector. 
Furthermore, the time derivative in x refers to the i-dependence of the points on this separatrix, which is determined 
by the deterministic equation of motion, i.e. Eq. (|^) with e = 0. Taking this and the definition (|^) of the current 
density into account, one obtains the general result 

T{t) ^-e J d'^-^n ■ DVp{x,t) . (11) 

5(t) 

As expected on naive grounds, a crossing of the separatrix is possible via a diffusive process only. 



III. PATH INTEGRALS AND THE SADDLE-POINT APPROXIMATION 



Let us briefly summarize the path integral description of the dynamics generated by a stochastic differential equa- 
tion (§). It is well-known from the literature (2|, |2|, |24[ ||, |6[ |^ ||, |9[ |^ |l|, ||, || that one can represent the 



conditional probability density p(xf , t/|xo, to) to find the particle at time tf at position xj if it has started at time to 
at position xq as weighted sum over all paths between these two points, 

x(«/)=x/ 

p(x^,i/|xo,to)= I Px(t)e-^[-Wl/^ . (12) 

x(to)=Xo 

Here, S'[x(t)] is the effective action or the Onsager-Machlup functional, i.e. 



S[x{t)] / dtL(x(i),i(i),0 , (13) 

with the corresponding Lagrangian 



to 



i(x, X, t) ■=\[^- F(x, t)] ■ D-i [x - F(x, t)] . (14) 
We remark that a prepoint-discretization scheme E3, Esj (not to be confused with the Ito-scheme in stochastic 



calculus [Q) has been implicitly adopted in the path integral (^2|) implying the measure det (47r e At D) . Other 
"discretization schemes" Q would give rise to a somewhat modified path-integral formalism but would — of 

course — lead to identical results as far as the actual stochastic dynamics (^ is concerned. 
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It should be clear that in general an exact evaluation of the path integral ( p^ is impossible. Nevertheless, it repre- 
sents an advantageous starting-point for a systematic weak-noise approximation: For small e, the main contributions 
to the path integral ( |l^ ) stem from the surrounding of action minimizing paths. Carrying out a infinite-dimensional 
saddle-point approximation around these paths, one arrives at 

g-0fc(x/,t/)/e 

piKf,tf\Ko,to)=J2 -[l + 0{e)], (15) 

where the sum runs over all paths x^, (i) that (locally) minimize the action (p^. These obey the Euler-Lagrange 
equations 

i^Ut) + [D {VF{^l{t),t)f D-i - VF(xI,(t),t)] ±Ut) - D {VF{^l{t),t)f Q-^ F(x^(t),i) - ^F(x*(t),<) ^ (16) 
together with the boundary conditions 

x^(to)=xo , ^litf)=^f. (17) 

Here and in the following, VF(x, i) denotes the Jacobian matrix with components (VF(x, := dFi{x,t)/dxj and 
the transpose of a matrix M is written as M"^. Introducing the canonical momentum 

f)T 1 

p:=^ = 2D-Mx-F(x,i)] , (18) 

we can, via the usual Legendre transformation, pass to the equivalent Hamiltonian dynamics, defined by the Hamil- 
tonian 

i?(x, p,i):=p-x-L = p- Dp + p- F(x, t) (19) 
and the corresponding canonical equations 

pUt) = -{'^Fi<it),t)fpl{t) ■ (20) 
The value of the action (j)k{xf,tf) of a path x^(i) is then given by 

M^f,tf) S[^Utf)] = rdtplit) ■ Dplit) , (21) 

J to 

where we have suppressed in favor of notational brevity the dependence on the initial condition xq at time Iq. It 
is noteworthy to point out that in Eq. (^9|) the momentum enters both quadratically and linearly: the latter linear 
contribution mimics a magnetic field-like, time-dependent vector potential contribution. For later use, we also recall 
the well-known result from classical mechanics [ [l4[ that the derivative of the extremal action with respect to its 
endpoint xj equals the canonical conjugate momentum at time t/, 

Plitf) = . (22) 



The yet unspecified prefactor term Q^(t/) in ( |15[ ) is given by the solution at time t/ of the following second order 
differential equation 

QUt) - VF(xI,(t), t) QUt) - QUt) QUt)-' D (VF(xI,(i),t))^ Q-^ QUt) 

VF(x^(t), t) D (VF(x^(i), t))^ - D (VF(x^ (t), t))^ D'' VF{^Ut),t) 

-^[VF(x^(i),i)-l-D(VF(x^(t),i))^D-i] +2(DpI,(t)).(VVF(x^(t),t))|Q^(0 = (23) 

with initial conditions 

QUto)^0, QUto)^D. (24) 
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Here, we have introduced the tensor of third order (WF(x, t))iji := 
with a matrix appearing in Eq. (123) is defined as 



d'^ Fi{x,t) /dxjdxi, whereby the scalar product 



[{Dplit}) ■ (VVF(x*(t),i))],, = Y.^Dpl{t)] 



dxidxj 



(25) 



The derivation of Eqs. ( p3| , p4[ ) in the framework of the time-discretized version of the path integral (|T^) proceeds along 
the same general lines as in the one-dimensional case presented in Ref. ]3^ . However, the rather technical details are 
beyond the scope of the present paper. 

Though in principle, the prefactor Q^(i/) is completely determined by Eqs. ( |23| , p^ ), it will turn out to be more 
convenient to consider the following quantity: 



D-i QUt) QWr' - {VF{Kl{t),t)f D-i - D-i VF(xI,(i), t) 



With the help of Eq. ( p3| ) one can then verify that G*f.{t) fulfills for t > to the matrix Riccati equation 

GUt) = ~ 2 GUt) D GUt) (VF(x^(i), t)f Gl{t) 

- GUt)VF{^l{t),t)~pl{t)-VVFi^Ut),t) . 



(26) 



(27) 



However, as can be inferred from Eqs. ( p^ , p6| ), there is no well defined initial condition for Gl{t) at time to. Instead, 
we have to content ourselves with the relations 



]im[Gl{t)Ql{t)] = -l and Gl{to)-' = 



By tracing over both sides of Eq. 
from Eq. (15): 



dt 



(28) 



Eq) one obtains a linear first order differential equation for the determinant det Q^. (t) 



detQlit) = 2Tr [(VF(x* (<), i))^ + DG^(t)] detQ*(i) . 



(29) 



The initial condition detQ^.(io) — follows immediately from Eq. (p4|). 

We conclude this section with two remarks on the above mentioned two different approaches for the calculation of 
the prefactor in Eq. (|l5|). First, we want to point out that one can identify G^(t/) with the Hessian of the action 
function 0^ (x^ , t j ) : 



Gfc(i/) 



9x/ 9x/ 



(30) 



The proof of this relation proceeds analogously to that presented in Ref. for the one-dimensional case, though 
the calculational details are by far more involved. Thus, Eq. (|2^) is equivalent to the matrix Riccati equation for the 
second derivatives of the action used elsewhere in the literature [Ml. However, while there this equation is derived 
by inserting a WKB-type ansatz in the Fokker-Planck equation (|), we entirely work here within the path integral 
formalism. In particular, we avoid problems commonly encountered due to the non-analytic nature of the WKB-action 
near the separatrix As a second remark, we note that an advantage of Eqs. (0, ^9|) over Eq. (|2^) lies in 

the fact that in the former set of equations no terms appear. Since later on we will be interested in the case of 
singular D, this fact presents a favorable feature from a technical point of view. 



IV. EVALUATION OF THE ESCAPE RATE FOR WEAK NOISE 



For small noise strengths e, we can now insert the saddle-point approximation (y_5|) into the rate expression (|l 
Using Eq. (^), this yields 

/■ ^ , _ P~<t>k{xl(tf).tf)/e 

k [det(4 7reQ*(t/))]^ 

Note that for the integrand the boundary condition ( p^ ) with Xf being the integration variable on the separatrix S{tf) 
is implicitly understood. Parametrizing the surface integral (|3l|) by the d — 1 dimensional vector s = (si, . . . , s^-i)^. 
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i.e. S{tf) =: {xsop(s, t/)|s € ^} for a certain subset -'^ C R'^ ^ of the d ~ 1 dimensional parameter space, we 
obtain 

where we have introduced the measure 

g(s, i/) := det -— . (33) 



u.u' — l....d—l 



For weak noise, the main contribution to the integral ( ^ ) comes again from the surrounding of the minima of the 
actions (f>ki^l.(tf),tf), and we can evaluate the integral in a saddle-point approximation. The corresponding extremal 
condition which determines the values then assumes the form, cf. also Eq. (22), 

dM-s.Asl,tf),tf) ^p.(,^). y,t,) (..!,. ..,d-l). (34) 



Consequently, the momentum p'^{tf) has to be perpendicular to the separatrix S{tf) at the end point Xscp(s^., tf). 
Furthermore, we need the Hessian of the action as a function of the parameter vector s. Via Eqs. ( p2| ) and (^6|) it can 
be expressed as 

scp(Sj,, 

«„* f)„* «„* k\ f) ^ * 

S^x (s* tA ^ ' 

+ Pfc(^/)- (z.,z/ =l,...,d-l) . 

Later on, we shall evaluate the last equation within a linearization around the unstable periodic orbit of our deter- 
ministic dynamics. In that case, the second term on the r.h.s of Eq. (35) vanishes. The determinant of the Hessian 
is thus given by the first term, which can be written as (cf. Appendix ): 



det { ^^^SSMM . GUtf) ^^£££MlM^ = 3(4, tj) nisltj) ■ Gl{tf)-' n(4, tf) det Gl{tf) 



(36) 



L'.L' 



Using the fact that the momentum pl{tf) is parallel to the normal vector n(sj!,t/) (cf. Eq. (|34|)), we then obtain for 
the rate (^2|) the important intermediate result 

r{tf)^^^=y , Pk[tf) uPk(tf)e [1 + 0(6)]. (37) 

V2^e^ ^pUtf)-Gl{tf)-^pUtf) detQlitf) detG*(t^) 



A. The action minimizing paths 

To proceed further in the evaluation of the rate formula (|37|), it is necessary to gain more insight about the nature 
of the action minimizing paths x^(t), whose dynamics is governed by the Hamiltonian equations ( pO| ) supplemented 
by the boundary conditions (|l7|). In view of our choice for the initial conditions (cf. discussion after Eq. (g|)) and our 
result (^7|), the latter assume the form 

x^.(to) = x,(io) , x^(i/) Xscp(Sfe,i/) , (38) 

where the parameter values are restricted by the relations (|3^). 

For given values of io and tf, in the generic case only one of the solutions x^(i) of this boundary problem will 
represent a global minimum of the action (|2l|). Denoting this path for the moment by x^ (t), it is clear that for large 
time differences tf — Iq owing to the form of the action ( p^ ) and ( p^ this path x| (t) will spent most of its time near a 
deterministic trajectory, i.e., x^(t) « F(x^(t),t). With regard to the boundary conditions (|38|), this means that after 
its start at time tg the path will closely follow the stable periodic orbit Xs{t) for some time. Subsequently it switches 
over into the vicinity of the separatrix S{t), thereby accumulating the main part of its action, to remain there until 
its end at time tf. With respect to the position of the end point on the separatrix, if we assume that ultimately all 
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deterministic trajectories on the separatrix converge to the unstable periodic orbit Xtj(t), the same will also hold true 
for the path (t) , which will thus end nearer and nearer to x„ (t) the more time it is able to spend in a close vicinity 
of the separatrix. Since the duration of the sojourns near the periodic orbits is long compared to that of the transition 
in between, the path x? (t) is often called an "instanton" in the literature. 

If we consider now the limiting case Iq ^ —oo and — *■ oo (in the following abbreviated as tf — to oo), there 
exists a well defined limit of x^ (t), in the sense that this path follows ever closer the periodic orbits, while retaining the 
shape of the intermediate segment. At the same time, its action Sfx^ ^t)] converges from above to a finite value. We 
observe that in the limit tf — to oo owing to the time-periodicity of the force field F(x, t) the action 5[x|(/: + nT)] is 
the same for all integers n. In other words, the global minimum of the action becomes countable infinitely degenerate. 
However, one can still safely assume that these minima are well separated in the space of all the paths appearing 
in (jl^), provided that the driving is neither too weak, too slow nor too fast. These limiting cases are thus not covered 
by our theory. We remind the reader of the situation in the static case, where the global minimum is also degenerate. 
However, while the degeneracy there is continuous in time (Goldstone mode), we are dealing here with a discrete 
degeneracy. 

As a consequence of the fact that the minimizing paths x^ (t) remain well separated, our rate- formula ( ^7|) becomes 
asymptotically exact for any (arbitrary but fixed) finite values of the driving amplitude and period as the noise strength 
e tends to zero. Apart from this fact that in the limit e — s- the 0{e) correction in the saddle point approximation 
( p^ and thus in ( ^7| ) vanishes, a more detailed quantitative statement seems difficult. On the other hand, for a given 
(small) noise strength e, we have to exclude extremely small driving amplitudes and extremely long or short driving 
periods since this would lead us effectively back to the static (undriven) escape problem, which requires a completely 
different treatment (especially of the (quasi-) Goldstone mode p7|, |3^ , ^ p5| , p6| ^ ) than in (|T5|). Put differently, 
in any of these three asymptotic regimes, the error 0{e) from ( |15| , |37| ) becomes very large. 

For the following, we introduce the symbol x*pj(i) for the limit of the path x^.(t) for tf — to — > oo, always keeping in 
mind that this path is only defined modulo time shifts by integer multiples of the driving period T. The corresponding 
action is defined analogously by 

0opt SiKpM ■ (39) 

Furthermore, all other quantities related to x.*p^{t) inherit the subscript "opt", for instance p*p(.(i) to name only one. 

Coming back to the case of finite times to and i/, we expect that as precursors of the limit tf — to ~* oo there 
exist besides the global minimum x? (i) additional, relative minima x^(t) with an only slightly larger action. After a 
suitable relabeling, each of them closely resembles an appropriately shifted "master path" x*pj (t + kT) (see Fig. |^) . 
For finite tf — to, we have a finite number of the order tf — tp/T of such paths. Thus, again without restriction of 
generality, we can assume that the sum in the rate formula ( p7|) runs from to a maximal value K{tf,to)' 

0<k<Kitf,to) = 0{{tf-to)/T) . (40) 

Especially, the path Xq (t) is the one that stays as long as possible in the vicinity of the stable periodic orbit x^ (t) and 
starts with its transition towards the separatrix at the latest possible moment. 



B. Linearization scheme around periodic orbits 



The discussion in the last section has shown that the paths x^,(i) spend most of their time near the periodic 
orbits x<,,u(t) of the deterministic dynamics. Thus, we may gain further insight in their behavior if we linearize the 
time-dependent force field around these orbits: 

F(x, t) « F(x,^„(t), t) + VF(x,,„(i), t)(x - x,,„(t)) . (41) 

Within this approximation, the Hamiltonian equations (^0|) assume the form 

A±l{t) = VF(x,,„(t), t) Ax^(i) + 2 D plit) 

where we have introduced the deviations from the periodic orbits, 

Ax*(t):=x*(t)-x,,„(0 . (43) 

As we shall see later, for the evaluation of the rate formula (^^, one only needs the solution for the momentum 
equation in the second line of Eq. (p^). According to Floquet's theory [Esl, this solution can be obtained from the 
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FIG. 3: Action minimizing paths xl{t), k — 0, . . . , K{tf,to) = 3 (solid) and corresponding master paths x*pt{t + kT) (dashed) 
between stable and unstable periodic orbits (dotted). Depicted is the one- dimensional, overdamped case d = 1 with an 
additively, harmonically driven piecewise parabolic potential from Ref. [pll p2| . 



eigenvalues and the left eigenvectors of the eigenvalue problem 



- + vr(x,,„(i),i) 



Kuit)^KuKuit) , *?.Ji + T) = *:,Ji) ,a = l,...,d 



(44) 



together with the corresponding right eigenvectors ^^'"(f), i.e. the solutions of Eq. ( p^ ) with (VF(xs,u(i), t))'^ . Upon 
proper normalization of these eigenvectors at equal times, i.e. 



we are able to write the solution for the momentum equation as 



(45) 



(46) 



where ti are arbitrary reference times for which the linearization (^ij) is valid. Since for the master path we have 
Ax*p^(i) ^ for i ^ ±oo, it follows with Eq. ( ^2| ) that limt^±oo p*pj(t) = 0. Consequently, we must require 
ReA" > (Re A" < 0) for all a with j,(ti) • p*pt(ii) ^ 0. In the following, we denote sums over this subset of 
eigenvectors by an apostrophe on the sum sign. 

With respect to the prefactor quantities G^(t) and detQ^(t), we first observe that within the linearization ( |4l| ) the 
matrix Riccati equation (P7|) assumes the form 



GUt) = -2 GUt) D Gl(t) - (VF(x,,„(t), t)f Gl{t) ~ GUt) VF(x,,„(i), t) . 
A transformation to the inverse matrix G^(<)^^ yields the linear matrix differential equation 

^ [Glit)-^] = 2 D + GI.(t)-i (VF(x,,„(i), t)f + VF(x,,„(t), <) Gl{t)-^ , 
where the initial condition is given by Eq. (|2^) as 

Gl{to)-'^0. 

Furthermore, multiplication of Eq. (^7|) with G^(i)^^ shows that 



Tr 



Gimiit)-' = -2Tr [(VF(x,,„(i),i))^ + D Gl{t)] 



(47) 
(48) 
(49) 
(50) 



With the help of the linearized version of Eq. (|29|) we then find for the last two terms of the product appearing in the 
denominator of our rate expression (B^ the result 



det Gl{t) det Ql{t) = const. —: fig^u 



(51) 
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Note that while is fixed for all paths ^%{t) by the initial condition (28) to the value fig = 2~ , /x„ depends both 
on the index k and on the explicit form of the time-dependent force field F(x, t). Later on, we only need /i„ for the 
master path, which we denote by 

//opt := lim det Q^^S) det Ql^,{t) ■ (52) 



The rest of the product in the denominator of (|37| ) can be rewritten by using the equation of motion ( |48| ) for G^(<) ^ 
together with the Hamiltonian equation ( ^ ) for \)%{t). This yields 

^[p^(t).G^(i)-ip^(t)] =2p;:(i).Dp:(i) (53) 

or equivalently in integrated form 

Pllt) ■ GUt)-' Plit) = plih) ■ GUh)-' pUh) + 2 fdt' pUt') ■ D plit') . (54) 

For the master path x*pj (t) , the integral on the right hand side of the last equation has to converge in the limit t ~* oo 
in order that the action (/)opt assumes a finite value. Hence, the quantity on the left hand side is also well defined in 
this limit and we obtain for t oo after a subsequent renaming tu t 

p:^,{t) ■ G:^,{t)~^ p:^,{t) = q,^, - p:^,{t) ■ A„(i) p:^,{t) (55) 

Here we have exploited Eq. ( ^ ) for the master path and furthermore introduced the abbreviations 

^opt ^hm p*pt(i) • Glp^{t)-^ Popt(^) : (56) 

, poo 

A„(i) := 2^ *^(t)*(^(t)^ / di'e-(^S+A;^)(*'-t)$t-(t') • D*t./3(i') . (57) 

Note that owing to the time-periodicity of the Floquet solutions, A„(t) is also a periodic function of t: 

Au{t + kT) = A„(i) for all times t. (58) 

C. Approximation in terms of the master path 

After having gained sufficient insight into the nature of the paths x^(t), we now come back to the evaluation of 
the rate expression (|37|). The main idea is to approxi mate a ll quantities related to the paths x^(t) in terms of the 



corresponding master path, since, as discussed in Sect. IVA, both resemble each other closely. However, while the 



boundary conditions for the master path are x*p^(t) — Xs(i) — > for t — > —00 and x*pj(t) — x„(i) — > for f — > 00, the 
paths x^ {t) have to satisfy Eq. (^8|) . From the results of our previous work , we know that while we can safely 
neglect these deviations at the initial time t^, the boundary at time tf has to be treated more carefully. First, we 
modify the upper boundary condition for x*p^(t) by requiring that 

tk:=tf + kT, ^k--^Kptitk) (59) 

are the new end time and end point, respectively. In other words, we simply truncate the path x*^^{t + kT) corre- 
sponding to x^(t) at the time tf. Obviously, this new path still satisfies the Hamiltonian equations (pG) and is thus 
an extremizing path. The value of its action follows from the definitions (^l|) and (|39| ) as 



dt p^pt (t) • D p^pt (i) = 0opt - dt pjpt {t) ■ D p:pt (t) . (60) 

-00 J tk 

Next, we express the action ( pl| ) of the path x^(t) by expanding the one belonging to the associated master path 
x*pt(t-|-fcT) in powers of the difference 6xl{tf) := x^(t^) — x*pt(tfc) = Ax^(iy) — Ax*pj(ifc) between their endpoints, 

M^scAsltj),tj) = </.opt(Xfc,t,.) + ^'^°P^(^^-'^^) . S^Utj) + ■■■ . (61) 
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The second term on the right hand side can be rewritten as p*pj(<fc) • Ax.l{tf) — p* ^(tfe) • Ax*pj(tfe). Here, the first 
scalar product is zero, as can been seen as follows: We know that the endpoint x^(i^) lies one the separatrix S{tf) 
and thus a deterministic solution Xdet(i) starting at time tk from Xdet(i/c) = will for all times stay on S{t). Its 

deviation from the unstable periodic orbit is denoted by Axdct (t) ■ Within our linearization ( ^l|) one can then readily 
verify with the help of the dynamical equation ( ^2| ) that d[p*p^(t) • Ax.dct{t)]/dt = 0. Hence, this scalar product is a 
constant of motion and we can infer that 

Pl^.itk) ■ A^litf) = p:^,{tk) ■ Axdct(tfe) = /nn [p:^,{t) ■ Axdet(0] = . (62) 



Here, the last equality follows from the boundedness of Xdct(0 near the saddle point (cf. the discussion below Eq. 
together with p*pj(t) ioi t ^ oo. 

The second scalar product, p*pj(tfc) • Ax*pj(tfc), can be determined by a similar consideration: The linearized 
Hamiltonian dynamics (^) yields 

^ [P^pt (t) ■ Ax:pt it)] =2 p^pt it) ■ D p:pt (t) , (63) 
and since again this scalar product has to vanish in the limit ^ — > oo, we obtain 

Pl^.ih) ■ Ax;pt(i,) = -2 / dt p;^,{t) ■ Dpl^,{t) . (64) 
Altogether, the approximation ( |6l| ) for the value of the action of the path x^(t) thus takes the form 

poo 

0fc(x,ep(s^i/),i/)=0opt+ / dtp:^,{t + kT)-Dp:^,{t + kT) + --- . (65) 

The fc-dependence of the exponential term in the rate expression ( p7|) is thus reduced to a simple sum of exponential 
factors as given by the time-dependence (^6|) of the integrand in (|65|). With respect to the prefactor terms in the 



rate expression (37), we already now from the results of our previous work |32[ | that up to corrections of the order 
^(|Popt(*fc)P) we can directly substitute all quantities belonging to the path x^(i) by those of the corresponding 
master path x*pj (t + kT) . Consistently, we also neglect the terms on the right hand side of Eq. (^^ , which are of the 
same order. 



D. Rate formula 



After the discussion in the previous section, it is just a matter of collecting everything in order to arrive at the 
central result of the present work, namely the asymptotic value of the instantaneous escape rate. 



T{t) ~ \/iaopte-^°p'/^ACopt(i,e) 



Here, we have introduced the quantities 



^opt 



[2''+\T^ \Ympl^,{t) ■ G:p,(t)-ip:pt(i) detQ:pt(i) detG:pt(0 

K{t,to 



u ^ ^ p*pt(t + /cr).Dp*pt(i + A;r) 

Kopt(t,£) ■.= T 2^ 



exp 



k=0 



dt' p:^,it' + kT)-Dp:^,it' + kT) 



(66) 

(67) 
(68) 



With the help of Eqs. ( ^ , |57| ) and the time-periodicity of the Floquet solutions (t) we can rewrite the last expression 

as 

K{t,t„} 



u \ -r \- bopt(t) • Bkity D Bk{t) bopt(t) 
Kopt(i, — ^ ^ cxp 

B,{t) := ^'e-^"'^^*t,a(i)$a(i)T ^ B,(t + r) , 

a 

bopt(i) := hm J2' *^ (i) • p:p,(£) *t'"(0 . 



^Opt 

{t) ■ Bkit)' K{t) Bfe(t)bopt(t) 
2e 



(69) 
(70) 
(71) 
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Since Bk(t) in ( |70| ) is a sum of terms that exponentially decrease with fc, there is a competition in the sum ( |69| ) 
between a pre-exponential factor which quickly decreases with k and an exponential term increasing with k. The 
main contribution to this sum thus comes from a few k values around a number k{t) that is for small noise strength e 
much larger than but at the same moment, for large enough tf — tg, still much smaller than K{tf, tp). Hence, up to 
an exponentially small error in e, we can extend the summation range in Eq. (|6^) and thus in Eq. ( |68[ ) to all integers 
k. This allows us to identify two important features of the time-dependence of Kopt(i, e)- To prove the first one, the 
time-periodicity 



■'opt 



{t + T,e) = Koptit,e) , (72) 



we merely have to shift the summation index in Eq. ( pq ) by 1. For asymptotic times and small noise strengths, the 
instantaneous escape rate ( |66| ) is thus, as expected, a periodic function of the time t. To establish the second property, 
we use furthermore that the pre-exponential term in Eq. (^8|) is just the negative time-derivative of the expression in 
the exponential. Again up to an exponentially small correction, we then obtain 



1 



t+T 

dt'Kopt(t',e) = 1 . (73) 



Inserting the asymptotic rate expression (|6^) into the definition (^) for the time-averaged rate and using the last 
identity, we thus arrive at the second main result of the present work, namely 

f ~ V^aopte-^-'/^ . (74) 

Hence, the noise-strength dependence of the time-averaged rate is of the form of an Arrhenius-type, exponentially 
leading term times an e-dependent prefactor. Note that for a system in thermal equilibrium i.e. for a time- independent 
force field F(x), the escape rate is given by an exponentially leading Arrhenius factor, which contains the barrier against 
the escape, multiplied by an e-independent prefactor, which depends only on local properties of the force field at the 
barrier and in the well |19| . In comparison with this equilibrium rate structure, we observe two crucial differences in 
([z^). First, in the present, non-equilibrium situation, both the effective potential barrier in the exponential and the 
prefactor depend in a non-trivial way on global properties of the force field F(x,t), and can thus in general only be 
determined by means of a numerical computation. Second, the non-equilibrium case exhibits an e-dependence of the 
prefactor. 



V. EXPLICIT RESULTS FOR A METASTABLE PIECEWISE PARABOLIC POTENTIAL 



While, in general, one has to resort to numerical methods for the evaluation of the rate expression (|66|), in the special 
case of a particle moving in a one-dimensional piecewise parabolic potential with an additive sinusoidal driving (^ a 
complete analytical treatment is possible. In the following, we work out the simplest such example with two parabolic 
pieces, and compare the so obtained analytical predictions with numerical results from a Monte-Carlo simulation of 
the stochastic dynamics (^). 

Let us thus consider the force field deriving from the piecewise parabolic potential of the form 

V(x<Q)^\mujl[{x-x,f -x^;\ 

\ (75) 
V{x>Q)=.-mujl[{x-Xuf-xl] , 

where Xs <Q and a;„ > denote the position of the potential minimum and maximum, respectively, with corresponding 
curvatures 

muol>Q and mujl<Q. (76) 

Note that within our notation uju is an imaginary number. To ensure the continuity of the corresponding force field 
F{x,t), 



F{x < 0,t) = -m,ujl [x -Xs)+ A sin(ili) 
F{x > 0, t) = -m Lol {x - Xu) + A sm{U. t) , 



- (77) 
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at the point a; = 0, we have to impose the additional restriction cu'^Xg = Selecting as independent model 

parameters lu'^, uj'^, and the static potential barrier AV :— V{xu) — V{xs), the fixed points Xg^u can be expressed 
through 



UJ?, x„ 



l2AV w?!^ 



y m ^2 + |a;„p 

If we require furthermore that the periodic orbits Xs.„(t) do not cross the point a; = 0, i.e. 

Xs{t) < < Xu{t) 

for all times t, which is granted if and only if the conditions 

A' < [7^ + in' - < J^] xl^ 

are fulfilled for both the "s" and the "u" indices, we obtain 

A 1^ cos{n t) + {n"^ -ujU sm{nt) 

Xs,u(t) = Xg II ■ — — ; - — . 

m 7^ 52^ + — 

For the determination of the master path x*^^ (t) , we use the specific form (^) of the force field and of the diffusion 
matrix to express the Hamiltonian equations of motion (Efl) in the form of the second order differential equations 



(78) 

(79) 
(80) 

(81) 



P*v,optit) - 7 P*v, optit) + ^luP*v, optit) = . 



(82) 



Here, the index "s" ("u") applies for all times t where ^pt(t) < (2:*pt(i) > 0). Furthermore, we have already set 
the strength S of the auxiliary noise term appearing in (|p to 0, since Eq. (|2|) and all the following relations remain 
well-defined for a singular diffusion matrix D. The system (^2|) of linear equations can be readily solved, if we restrict 
ourselves to the case where the master path xlp^{t) crosses the point a; = exactly once, say at time ti, 



t = ti 



(83) 



Notice that owing to the form of the Hamiltonian equations ( p2|) both x*p^ {t) and p.* (t) have to be continuously 
differentiable at this time ti. For the following considerations it is convenient to introduce the frequencies 



S.U 1 



(84) 



which are the Floquet eigenvalues from Eq. (Q). We observe that A,j < 0, such that the respective terms do not 
appear in the sums (p7|j70|j7l|) . The corresponding Floquet states are given by 



As u 



and *t.±(i) 



+7 

1 



Hence, the matrices Bfc(i) from Eq. ( |70| ) take the form 



Bfe(t) = 



-,-XZkT 



K+1 



(85) 



(86) 



Taking into account the boundary conditions Ax*p^{t) — s- and p* opt('') ~^ for t ±oo, we obtain for the solutions 
ofEq. 



A<pt(i <ti) 
A<pt(i>ti) 



1 



As — Ac 



^pJ;,opt(^i) - A^ xsih) 



-AJ-(t-ti) 



— K,opt(^l) - As Xsih) 



-A+-(t-ti) 



1 



1 



P*v,optih) - A+ Xuih) 



m ^ J 



,opt 



TO 



A? - A 



(87) 
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Next one can infer from the above mentioned requirement of a continuous first derivative at time t of both x*^^ (t) 
and pl^opti^) the two relations 



K,opt(*i) ^rn^X+Xsih) 



(88) 



These relations can be used to fix the two so far unknown quantities ti and p* opt(^i) the above expressions: A 
straightforward, but somewhat tedious calculation yields 



UJg 



cos(r2 ti) 



(89) 
(90) 



Obviously, Eq. ( |89| ) has two solutions within each driving period T. We anticipate that only one corresponds to a 
minimum of the action, and hence to the master path. Thus, we fix ti (up to the usual degeneracy under t t + T) 
by requiring additionally that 



A cos(riii) 



> 



Combining (^-(|9l]) we arrive at 



* /X > -I- 1^1 



Xu 



where we have introduced the frequency j/ via 



x:, 



1/2 



(91) 



(92) 



(93) 



To conclude the discussion of the master path, we note that in order to check the consistency of the solution ( p7| ) with 
the condition (^) it is necessary to solve a transcendental equation, which has to be done numerically. Qualitatively, 
one expects a break-down of Eq. ( p3[ ) in the deterministically underdamped regime, i.e. 7/2 < Wg, where the path 
{xopt{t), Vopt{t)) in phase space leaves the vicinity of the stable periodic orbit by way of a spiraling orbit which crosses 
the line a; = several times. 



Inserting p* ^ ^ {t) from Eq. (p^ and 



into the definition (tel^,p9) yields the action of the master path: 



"opt 



AV 



1/2' 



(94) 



The periodic driving thus leads to an "effective potential barrier" (popt which is smaller than the static barrier AV 
to which it reduces in the limits A ^ or SI — > c». While (popt is monotonically decreasing with increasing driving 
amplitude A, the dependency on the driving frequency Q is more complicated. In particular one observes resonance 
behavior near the frequency ojg of the bottom well. Had we chosen the opposite inequality in Eq. (|9l| ) , a plus instead 
of the minus sign would have appeared in the Eq. (|94|). Thus, the condition (91) does indeed single out the minimum 
of the action and thus the desired solution for the master path x*^ (t) . Using the result for pJJ^opt(*i) from Eq. (|9^), 
we obtain the following explicit expression for bopt(i) from Eq. (O): 



bopt(0 = 



1 



)K,opt(^i)e-^"(*-*^) 



(95) 



Finally, we turn to the determination of the prefactor quantity aopt ■ We first note that in the present case the linear 
differential equation ( ^ ) does not only hold approximatively but is exact for all times t ^ ti. Owing to to —00, 
the solution has approached its stationary value, and we obtain 



G:pt(i < h) 



raujj. 

TO 



(96) 
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and hence with Eq. 



det Qlpdt<h) 



1 



(97) 



At this point, we make use of the fact that detQ*pt(i) is continuous at time ti according to Eq. (|29|), while G*pt(t) 
jumps due to Eq. (p7|). With d'^F{x,t)/dx'^ = m (w^ — S{x) one thus can infer that 



, I 12\ ™^;pt(*i)-P:.opt(*i) I 12 n' 
+ \^u\ ) ^T^) \cJu\ 

m 



(98) 



This yields together with Eq. (51), which again is exact for the piecewise parabolic potential, the intermediate result 



^hmdetG:p,(i) detQl^,{t) = 



f 2 , \ \2\ ™^opt(^l) Pv.opti^i) I |2 

[uj^ + \uju\ ) — - — , m\uj, 



(99) 



Applying two times Eq. (jsj), one time for t ^ to ^ —oo using the initial condition (^) and a second time for t ^ oo, 
yields after some algebra another necessary quantity, namely 



hm p:pt(i)-G:pt(0"'p:pt(i) 



,2,1,, |2\ ™iopt(^l) -P^,opt(^l) I ,2' 



2 (j)o 



pt 



(100) 



Inserting these expressions into the definition ( |67|) of Q;opt and using Eqs. (|2C|), ( |92| ) and { p4\) we arrive at 



-*-opt 



^^^2 '^^opt(^l) -P^,opt(^l) 



TOX* t(tl) 



•'opt 



-1/2 



(101) 



Again, our choice for the sign in Eq. (91) is justified a posteriori, since it guarantees that m x*p^(ti) — p* opt(^i) is a 
positive quantity. Eventually, together with Eqs. ([78|), (pl|), (p9|), ( pO[ ) and (p^), the last expression can be rewritten 
as 



a, 



opt 



1^1 



2mAy jy** 



7^ 



pt 



1/2 



(102) 



We remark that the structure of the results ( p4[ ) and (102) closely resembles those in the over damp ed case, as obtained 
in Ref. js^, |3|]. In particular, this limiting case m ^ is correctly reproduced by (|^ and ( |l02| ). 

We have checked our above analytical predictions by comparing them with results from Monte-Carlo simulations 
of the stochastic dynamics (^), which yield the mean time iexit necessary for an exit out of the driven well. The 
time-averaged rate is related to this quantity via |49] 



(103) 



Since the necessary simulation times diverge exponentially with decreasing noise strength e (cf. ( |74|) and (lOS)), the 
numerical determination of the mean exit time becomes impossible for extremely small e. For similar reasons, the 
numerical determination of the time-resolved rates T{t) is ruled out. 

Figures ^-|^ depict the time-averaged rate f as a function of various system parameters. Parameter regions for 
which condition ( ^3|) is not fulfilled are indicated, which in particular applies for the case of weak damping rj and the 



vicinity of the deterministic resonance at f2 



7 /2, as discussed above. Outside of these regions, already for a 



noise strength of e = 1, the agreement between theory and simulation is very good (compare also inset of Fig. 1^). 

Furthermore, as one can see in Fig. ^, the theory breaks down both for small and large driving frequencies if the 
noise strength e is kept fixed. However, for e — > 0, we again observe the convergence of numerically determined rate 
towards th e the oretical approximation (cf. inset of Fig. ^). This behavior is in full accordance with our predictions 
from Sect. lIV A|. 
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FIG. 4: Time-averaged rate f vs. driving frequency f2 for Xs — co^ — —1, Xu ~ tt;: 
(dimensionless units). Crosses: Inverse mean first exit time from simulations of Eq. 
mean value). Solid (dotted) line: Analytical prediction (j?^) (condition (js^) not fulfillec 
e — > for fixed driving frequency Q = 5 (all other parameters as in main panel). 



= m = ri = AV ^ A = l,e = 0.1 
(error bars: mean square deviation of 
. Inset: Convergence behavior off for 



1^ 10"* r 




FIG. 5: Time-averaged rate f vs. particle mass m for Xs = muj^ = —1, Xu = mtj? = rj = AV — A = Q, — l,e = 0.1 
(dimensionless units). Crosses: Inverse mean first exit time from simulations of Eq. (H) (error bars: mean square deviation 
of mean value). Solid (dotted) line: Analytical prediction ( [74|) ( condition ( |83| ) not fulfilled). Inset: Magnification of small m 
regime. Note that the deterministic resonance lies at m = \/l/2 ~ 0.707. 



VI. CONCLUSIONS 



With the present work, we have generalized earlier results for the noise-activated escape in time-periodically driven 
one-dimensional overdamped systems to an arbitrary number of dimensions. The basic idea of the path-integral 
approach put forward is that it is necessary to sum over all local, nearby minima of the relevant action that contribute 
to the escape rate. For asymptotically weak noise, these minima are well separated in the space of all paths, and the 
contribution of each of them can be obtained by a standard saddle-point approximation of the path-integral. This 
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FIG. 6: Time-averaged rate f vs. damping coefficient 7 = rj/m for Xa = (-t^t = — Ij ^u-r^ loI ~ m = W — A = Q, = l,e = 0.1 
(dimensionless units) . Crosses: Inverse mean first exit time from simulations of Eq. (p|) (error bars: mean square deviation of 
mean value). Solid (dotted) line: Analytical prediction (^) (condition ( |83| ) not fulfilled). Inset: Relative difference between 
analytical prediction and numerical result as a function of 7. 



behavior is in strong contrast to the undriven situation, where the occurrence of a (quasi-) Goldstone mode requires 
a more sophisticated treatment of the associated prefactor. 

Our central results express the time-instantaneous rate (|6^) and the time-averaged rate ( fz^ ) in terms of quantities 
belonging to a master path, which has in general to be obtained numerically as the solution of a minimization problem. 
However, the structure of the result already reveals two noteworthy, salient differences compared to the static case: 
both the exponentially leading Arrhenius factor and the prefactor not only depend on the relative barrier height and 
the local properties of the potential inside the well and at the barrier region, respectively, but in addition these two 
quantities depend as well sensitively on the global behavior of the time-dependent metastable potential. Additionally, 
a -^/e-dependence of the rate prefactor appears as a consequence of the non-equilibrium situation. Finally, we have 
been able to derive closed analytical rate expressions for a distinctive case, namely the sinusoidally driven Kramers 
problem with a metastable potential formed by two parabolas. 

We also like to point out here that the treatment of driven escape in meta-stable potential landscapes by use of 
path integral methods is close in spirit to prominent work by Jozef T. Devreese wherein he pioneered the challenging 
problem of polarons in magnetic fields |l|: in both cases the effective Hamiltonian involves a vector-potential 
like contribution which is linear in the canonical momentum variable, see Eq. (|l9|). Then, as pointed out repeatedly 
by Jozef, the standard variational principle based on the Feynman- Jensen inequality no longer applies ||5l|, but 
requires instead an appropriate modification |52 . 
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APPENDIX 



In this appendix, we shall prove Eq. (|3^) . Let us first abbreviate the d x d — 1 matrix S composed of the vectors 
9xsop(s^, t/)/9s^ ^, V — 1, ... ^d — 1 as columns. For notational brevity, we suppress here and in the following the 
asterix, all indices k, and the time arguments. Indicating furthermore with S\ « = 1, . . . , c?, the d — 1 x d — 1 matrix, 
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which consists of S without row i and with G*, i,j = 1, . . . , d, the matrix without row i and column j, we can apply 
a generalized version of the determinant multiplication theorem to the left hand side of Eq. ( ^6| ) , yielding 

det(S G S) = ^ det dot G] det S* . (A.l) 

Using a standard relation between the cofactors and the inverse of a matrix [S^ 

(-1)*+^' det G; = (G^i)j. det G , (A.2) 

we obtain the intermediate result 

det(SGS) =u-G"iudetG . (A.3) 

Here, we have introduced the rf-dimensional vector m := (— 1)*+'* det S\ i — 1, . . . , d. A comparison with Eq. (|3^ ) 
shows that the proof is complete, if we are able to demonstrate that u = y^n. To this end, we define another 
auxiliary quantity, namely the d x d matrix M := (S|n). Since the normal vector n of the separatrix S is orthogonal 
to all columns of S, the determinant of M is given by (cf. also Eq. (B3[)) 



det M = VdctM^M = Vdet S^S = ^5 (A-4) 
For the same reason, one has {M^^)di — Mid = nt, i = 1, . . . ,d, and thus 

= (-1)^+'^ det = det M iM-^)d^ = ^gn, . (A.5) 
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